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CURVE  FIT  TECHNIQUE  FOR  A  SMOOTH 
CURVE  USING  GAUSSIAN  SECTIONS 


INTRODUCTION 

In  optics.  It  Is  often  necessary  to  determine  the  total  luminous  flux  of  a  beam 
by  use  of  sampling  beam  receivers.  If  the  beam  Is  circular  In  cross  section  and 
It  Is  desirable  to  approximate  this  with  a  Gaussian  energy  distribution  along 
any  diameter,  then  It  Is  referred  to  as  Gaussian  distribution.  Receivers  placed 
on  a  diameter  of  this  beam  provide  discrete  samples  so  that  the  total  amount  of 
flux  may  be  calculated.  Thus,  a  technique  that  forms  a  smooth  Gaussian  patch 
between  any  number  of  discrete  sample  points  would  be  useful  for  a  total 
Incoming  flux  approximation.  This  report  deals  with  the  development  of  this 
technique. 

DISCUSSION 

Development  of  First  Gaussian  Section 

To  simplify  the  fitting  technique,  two  sampling  cases  are  assumed.  One  Is 
that  one  sample  point  has  an  Intensity  greater  than  all  other  sample  points,  and 
the  other  Is  that  two  adjacent  sample  points  have  the  same  luminous  Intensity. 
These  sampling  cases  lead  to  two  different  beginnings  for  the  Gaussian  sections 
surrounding  the  peak  point  and  will  be  dealt  with  separately. 

If  there  Is  a  single  peak  point,  represented  by  (xQ,  y0),  the  Gaussian  sec¬ 
tion  that  spans  to  the  next  point  on  either  side  must  be  found.  A  Gaussian  sec¬ 
tion  (Gn)  has  the  form 


-uxk  :.v 


Equation  (1) 


(xk»  yk)  «r«  the  coordinates  of  any  point  which  the  section 
passes  through 

un  Is  the  mean  point  of  the  Gaussian  section 
sn  Is  the  standard  deviation  of  the  Gaussian  section 
cn  is  the  normalizing  constant  for  the  Gaussian  section. 

The  peak  point  of  any  Gaussian  has  a  slope  of  0,  therefore,  the  first  section 
must  have  Its  peak  at  (Xg,  yg).  Here,  the  exponential  is  at  its  maximum,  1, 
therefore  the  exponent  is  0,  leading  us  to  realize  that 

Uj^  =  Xg  Equation  (2) 

Now,  since  the  exponential  has  a  value  of  1,  the  Gaussian  simplifies  to 

Cj  *  yQ  Equation  (3) 

Rearranging  (1)  for  the  standard  deviation  and  letting  (xi.yi)  be  the  other 
boundary  the  Gaussian  passes  through,  we  have 


Equation  (4) 


and  equations  (2)  through  (4)  describe  the  first  Gaussian  section. 

If  two  points  have  the  same  Intensity,  they  can  be  assumed  to  be  boundaries  of 
the  same  Gaussian.  Then  the  mean  of  the  section  is  the  midpoint  between  the  two 
sample  points.  Here  (x^,  y^)  represents  a  second  point  with  the  same  intensity 
as  the  first,  and  It  is  realized  that 


U1  *  (xi  +  xn)/2 


Equation  (5) 


To  simplify  matters,  an  assumption  will  be  made  that  the  measured  Intensities  at 
these  points  are  90%  of  the  peak  Intensity.  This  means  that 

C1  *  yj/0.9  «  y0/0.9  Equation  (6) 

By  using  equation  (4),  can  now  be  found.  Thus,  equations  (4)  through  (6) 
decribe  the  first  Gaussian  section  in  the  second  case. 

Development  of  Subsequent  Sections 

After  the  first  Gaussian  section  Is  fitted,  subsequent  Gaussian  sections 
must  be  fitted  to  assure  a  smooth  curve.  A  smooth  curve  Is  one  where  both  func¬ 
tions  have  equal  values  and  equal  slopes  at  the  point  of  junction.  This  con¬ 
dition  on  the  second  section  leads  to 

G'n(xk)  *  G'n+l(*k)  Equation  (7) 

Gn  Is  the  Gaussian  section  ending  at  (xfc,  yk)  and  Gn+l  Is  the  Gaussian  section 
extending  from  (x|<,  yk)  to  (xk+1,  yk+l).  Taking  the  derivative  of  Gn  from  (1) 
at  the  common  boundary  point  (xfc,  yfc)  we  have 


Gn'(xk) 


-yk 


(xk  -  un) 


Equation  (8) 


Performing  the  same  derivative  for  the  Gaussian  Gn+i  at  (xk,  yk)  we  have 


(xk  -  un+l) 


6n+l'(xk)  ■  -yk 


Equation  (9) 


Using  equation  (7),  equations  (8)  and  (9)  may  be  set  equal  to  each  other.  From 
this,  the  -yk  factors  may  he  cancelled  out  and  the  equation  may  then  be 
rearranged  to  solve  for  sn+i 


(«n+l)2 


(xk  -  un+i) 

(xk  -  un) 


Equation  (10) 


which  relates  sn+i  to  $n.  However,  this  equation  has  the  term  un+l  which  has 
yet  to  be  found.  To  do  so,  we  take  the  Gaussian  Gn+1  previously  defined  as 
being  bounded  by  the  points  (xk,  yk)  and  (xk+1,  yk+1).  The  above  Gaussian  must 
then  satisfy  the  following  condition: 


n+1  *yke 


Equation  (11) 


From  this  we  take  the  natural  logarithm  of  both  sides  of  the  equation  and 
rearrange  it  to  Isolate  un+i 


2(sn+i)  (In  yk+1  -  In  yk)  - 

(xk  "  Vl)2  *  (xk+l  "  %+i)2 


Equation  (12) 


Substituting  the  general  form  for  Sn+1  (from  equation  (10)),  and  rearranging 
terms  we  have 


2xkSn2(ln  yk+l  -  In  yk)  -  (xk2  -  xk+l2)(xk  -  un) 
2sn2Un  yk+1  -  In  yk)  -  2(xk  -  xk+l)(xk  -  un) 


Equation  (13) 


This  allows  un+i  to  be  found  solely  In  terms  of  known  quantities.  Then  un+l  may 
be  substituted  Into  equation  (10)  to  find  Sn+1*  These.  In  turn,  may  be  used  In 
equation  (1),  along  with  the  point  (x|<,  yi<),  to  solve  for  cn+i.  Thus  we  have  a 
method  of  finding  a  Gaussian  solely  from  the  preceding  Gaussian  and  its  boundary 
points.  This  method  will  work  for  any  case  and  Is  therefore  general. 

Computer  Testing  (Sample  run,  figures  1  and  2) 

To  facilitate  validation  of  the  curve-fit  technique  developed,  a  computer 
program  was  written  to  perform  the  task  of  finding  the  Gaussian  sections  of  a 
curve-fit  problem  by  the  methods  given  here.  (The  program,  run  on  a  Tektronix 
GS-4052  System,  is  listed  In  figure  2.) 

The  data  used  for  the  sample  run  consisted  of  the  peak,  90%,  50%,  and  10%  points 
of  a  Gaussian  curve  with  u  *  0,  s  «  1,  c  ■  0.398942.  The  printout  from  a  sample 
run  with  first  one  side's  curve  fit  and  then  the  other  side's  curve  fit  is  shown 
In  figure  1. 


Section 

Mean 

Stan.  Dev 

Norm  Cons 

1 

0 

0.999834481528 

0.398942 

2 

0 

0.999668990452 

0.398962878929 

3 

-4.415707426E-4 

0.999712877896 

0.399308591275 

Section 

Mean 

Stan.  Dev 

Norm  Cons 

1  0  0.999834481528 

2  0  0.999668990452 

3  4.415707426E-4  0.999712877896 


0.398942 

0.398962878929 

0.399308591275 


Figure  1.  Output  from  sample  curve-fit  program. 
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Figure  2.  Listing  of  sample  curve-fit  program. 
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The  sample  run  output  seems  quite  accurate.  To  determine  the  accuracy,  a  graph 
was  prepared  calculating  the  percent  deviation  of  the  curve-fitted  curve  with 
the  original  curve,  and  the  results  are  shown  in  figure  3. 


Figure  3.  Graph  of  percent  deviation  as  a  function  of  x. 

The  greatest  deviation  is  under  .15%,  an  error  small  enough  to  be  attributed  to 
computer  rounding  and  sample  error. 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  above  described  technique  has  led  to  a  deviation  from  a  norm  of  no  greater 
than  .15%.  This  devation  is  of  a  small  enough  magnitude  to  validate  the 
developed  method  of  curve-fitting.  Furthermore,  the  algorithm  that  does  the 
fitting  is  simple  enough  to  be  used  on  a  programmable  calculator. 


